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We use a standard bead-spring model and molecular dynamics simulations to study the static 
properties of symmetric linear multiblock copolymer chains and their blocks under poor solvent 
conditions in a dilute solution from the regime close to theta conditions, where the chains adopt 
a coil-like formation, to the poorer solvent regime where the chains collapse obtaining a globular 
formation and phase separation between the blocks occurs. We choose interaction parameters as 
is done for a standard model, i.e., the Lennard- Jones fluid and we consider symmetric chains, i.e., 
the multiblock copolymer consists of an even number n of alternating chemically different A and B 
blocks of the same length Na = Nb = N. We show how usual static properties of the individual 
blocks and the whole multiblock chain can reflect the phase behavior of such macromolecules. Also, 
how parameters, such as the number of blocks n can affect properties of the individual blocks, 
when chains are in a poor solvent for a certain range of n. A detailed discussion of the static 
properties of these symmetric multiblock copolymers is also given. Our results in combination 
with recent simulation results on the behavior of multiblock copolymer chains provide a complete 
picture for the behavior of these macromolecules under poor solvent conditions, at least for this 
most symmetrical case. Due to the standard choice of our parameters, our system can be used as a 
benchmark for related models, which aim at capturing the basic aspects of the behavior of various 
biological systems. 

PACS numbers: 87.15.A, 87.15. N, 87.15. B, 82.35.Pq, 82.35.Jk 
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I. INTRODUCTION 

Polymers combining chem.ically different mononicric 
units on their structure found very much theoretical and 
experimental interest, as, for instance, diblock copoly- 
mers, where in the case of a linear chain one block of the 
chain consists of one type of monomers, while the other 
part of another type of monomers. These macromolecules 
have already found particular interest as they form dif- 
ferent structures depending on their structural param- 
eters, thermodynamics conditions, etc. For a standard 
theoretical review one could see Ref. [l[ and references 
therein. We underline that the motivation for studying 
such systems stems from the fact that these materials 
are used in industrial scale [l| and also in very advanced 
technological applications (as an example one may look 
a very recent review Ref. [^). As an extension, regarding 
possible complexity in the structure or other aspects, of 
diblock copolymers is regarded multiblock copolymers, 
where in their simplest case more than two blocks com- 
prise a linear polymeric chain and the blocks are com- 
posed either of type A or type B monomers. The phase 
diagram of multiblock copolymer melts Q resembles the 
phase diagram of diblock copolymer chains [^4y| : where 
the scaling parameter xN controls the phase behavior of 
the system, as it has also been discussed in simulation 
studies H- For the most symmetrical case the lamel- 



lar structure is predicted, an approach that is particu- 
lar valid in the strong segregation limit. The case of 
multiblock copolymer chains is particularly interesting 
even when only two type of blocks (A, B) alternating 
along a linear chain in an infinite dilute solution is con- 
sidered |8l-[l5l|. In this case, it is sufficient to study the 
behavior of an isolated chain taking into account only 
interactions (often of short range) within the chain. A 
model describing a multiblock copolymer chain could also 
exhibit some interest due to its close relation to vari- 
ous toy- models |l6|, which try to mimic the behavior 
of various biomacromolecules in several processes [l7| . 
Such biomacromolecules are rather formed by periodi- 
cally repeated structural units ( "monomers" ) along their 
chain. Simple (coarse-grained) models jlSj incorporating 
the key aspects for a specific biomacromolecule are able 
to capture efficiently their relevant properties for which 
the model is designed. 

The phase behavior of a single multiblock copolymer 
chain has been extensively studied by means of molec- 
ular dynamics (MD) simulations by using a standard 
bead spring-model and considering interactions between 
beads according to a standard model (Lennard- Jones 
(LJ) fluid) J13l - [l5l |. Of course the use of MD simulations 
imposes certain limits in the range that the parameters 
may vary, for example, the range of the temperature and 
the chain length. However, it has been shown that the 



different regimes of interest for the phase behavior can ah 
be accessed by simulations [i-[lS [iS-El [11 ■ Multiblock 
copolymers composed of two different types of blocks (A 
and B) alternating along the linear chain and length of a 
block A being equal to the length of block B that adopt 
coil-like structures in agood solvent or at temperatures 
close to the Q \ViHl3i [l9| have been studied. The chem- 
ical difference of monomers should actually be kept re- 
sponsible for an expansion in the chain dimensions with 
respect to the equivalent homopolymer chains (same to- 
tal number of monomers) under the same thermodynamic 
conditions [ll|, ll2| . Below a certain temperature [ij] a 
multiblock copolymer chain collapses due to the incom- 
patibility between monomers and the solvent obtaining 
the so-called globule formation and different scenarios of 
phase behavior have been discussed [IJ, [1^ . Simulations 
have confirmed the validity of the parameter xN in defin- 
ing the phase behavior of a di/multi-block copolymer sys- 
tem, but it has also been discussed that the geometry of 
the microphase separated regions is controlled by the to- 
tal number of blocks n, as well as other parameters, i.e. 
the relative size and arrangement of the blocks [3] . The 
same has been shown to be true for a single multiblock 
copolymer chain under poor solvent conditions |13l4l5l |. 
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The phase behavior of such macromolecules has been 
discussed within the frame of the analysis of the formed 
clusters (which contain two or more blocks of the same 
type) and of the probability that a cluster of one type 
of blocks contains all A blocks while another cluster con- 
tains all B blocks J13l - [la |. i.e. a single lamellar struc- 
ture with an A-B interface between blocks of type A and 
blocks of type B. Thus, when this probability is unity 
we fall into the above single lamellar case. Then, for the 
range of temperatures (x oc Ae/T, Ae = const) and of 
chain lengths {nN, where N is the block length; blocks A 
and B have the same length) accessible to simulations this 
type of phase behavior has been seen for small number of 
blocks n and for values of block length A^ > 20. Another 
scenario suggests that the above type of phase behavior 
takes place with a certain probability {Pn^i=2 < 1, Nd 
being the number of clusters; such probability can be 
very high or very small depending on the values of N, n, 
and T). The third scenario corresponds to the case that 
the probability P/v^,=2 = (the formation of only two 
clusters with monomers of different type never occurs) 
and a symmetric variation in the number of clusters Nd 
around an average value (2 < N^i < n,n is the total num- 
ber of blocks of the chain) is seen. Of course, the analysis 
of the probability is performed separately for clusters of 
type A monomers and for clusters of type B monomers. 
Otherwise the A and B beads that belong to the same 
chain would always belong to a single cluster contain- 
ing all the beads of the chain. Similar conclusions have 
been also drawn by a similar model |8l-[l(j|. where also 
the geometry of these cluster phases has been analyzed 
in detail. Guided by the theory Q, one would rather 
expect in the long chain limit that a ground-state-type 
structure would be a single lamellar domain, where an in- 
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FIG. 1: (Color online) Definition of structural 
parameters describing our linear multiblock copolymer 
chains, n (in this case n = 4) is the number of different 
blocks A and B {ua = ns = 7^/2) denoted with different 
color (or grey tone) and N is the length of each block. 
All the blocks, irrespective of whether they are of type 
A or B, have the same length N. Then the total length 
of the chain is nN. 



terface between all A- and B-typc blocks is formed, simi- 
larly to what is known for multiblock copolymer melts for 
symmetric choice of parameters. Such a structure would 
have much less (unfavorable) A-B contacts compared to 
a multidomain structure of A and B clusters, which is 
kinctically favored in simulations. Recently, the static 
properties of the formed clusters for such systems have 
been discussed [IJl , where one is also able to detect in ad- 
dition to the phase boundaries the globule-coil crossover 
and the proximity of the O temperature through such an 
analysis. In this manuscript, we will show further, that 
properties such as the bond correlation and asphericity 
can provide a rather accurate evidence for detection of 
full phase separation (single lamellar structure) in multi- 
block copolymer chains. Additionally, we will argue that 
these properties can provide information on the length 
N and the number of blocks n that this phase behavior 
is taking place. 

Furthermore, we provide a complete understanding 
and a broad discussion of the static properties of such 
macromolecules under poor solvent conditions for the 
blocks and the multiblock chain as a whole, at least for 
the most symmetric case that is considered here, and 
for chain lengths and thermodynamic conditions accessi- 
ble to our simulations. As we have chosen in our study 
the same potential parameters as for the standard LJ 
fiuid [201, our results could also serve as a benchmark 
for similar polymeric and biological systems. We discuss 
our results only in the poor solvent regime, which is the 
most interesting case due to the different phase behavior 
that has been observed when the chain is in the collapsed 
state (globule state) |13l - [l5| . 

The rest of the paper is organized as follows. In Sec.HIl 
we describe our simulation model and the methods to 
analyze our results. In Sec. IIIII we present our results 
for the whole multiblock chain and the individual blocks 
and we discuss particular details of the studied proper- 
ties pertinent to the phase behavior of the multiblock 
chains. Finally, in Sec. II VI we summarize our conclusions 
in connection with recent results on the phase behavior 
of symmetric linear multiblock copolymers. 



II. MODEL AND METHODS TO ANALYZE 
THE SIMULATIONS 

As a detailed study of the systems under consideration 
is computationally demanding for reasons that will be 
discussed below, we have focused only on symmetric lin- 
ear multiblock copolymers, where blocks of type A beads 
and blocks of type B beads alternate along the chain. 
Then the length of all blocks, irrespective of whether be- 
ing of type A or B, is N, the number of A-type blocks 
UA is equal to the number of B-type blocks ub, with 
n = UA + TiB an even number denoting the total num- 
ber of blocks. A multiblock copolymer chain is described 
by the parameters shown schematically in Fig. [Jl where 
the different colors correspond to the different types of 
blocks, which alternate along the chain. 

To model our chains we use the standard bead-spring 
model [2l| - |26| to describe a multiblock copolymer chain, 
where all beads interact with a truncated and shifted 
LJ potential ULj{r) and nearest neighbors bonded to- 
gether along the chain also experience the finitely exten- 
sible nonlinear elastic potential UpENEif), ^ being the 
distance between the beads. Thus, 
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where a, j3 = A^B denote the different type of monomers, 
and the constant C is defined such that the potential is 
continuous at the cut-off {vc = 2.5). Henceforth units are 



chosen such that £° 



1, tr"" = 1, the Boltzmann con- 



stant fcs = 1, and in addition the mass w^j of beads 



is chosen to be unity. 
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For simplicity, cr^j ~ 1, but 



e>'V = eyy = leYfj = 1, in order to create an unmix- 
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ing tendency between monomers A and B belonging to 
different blocks as done in previous studies [ISl - llSl . p H - 
[29j . and as is used for a standard system (LJ fluid) J2G |. 
Therefore, Ae = eff - \l2{e^j + eff ) was kept the 
same throughout our simulations and x (x °^ ^^/T) was 
mainly varied by tuning the temperature T . The poten- 
tial of Eq. ([1]) acts between any pair of beads, irrespective 
of whether they are bonded or not. For bonded beads ad- 
ditionally the potential UFENEir) acts. 
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fcr^ln[l-(r/ro)2], < r < tq. (2) 



where UFENEif ^ fo) = oo, and hence tq is the maximal 
distance between bonded beads. We use the standard 
choice To = 1.5 and k = 30. We recall that for lin- 
ear chains the temperature for the present model has 
been roughly estimated as Tq si 3.0 (a more accurate 
value [30| Tq w 3.18 could only be established for chain 
length nN > 200, n = 1). Thus, we have considered 
in our study the range of temperatures 1.5 < T < 3.0. 
Moreover, in previous discussion of this system [l3| we 
showed that the chain enters in the globule regime at tem- 
perature T « 2.4. Therefore, phase separation between 
A and B monomers takes place below this temperature. 



For homopolymer chains this boundary is sharper com- 
pared to the case of multiblock copolymers where also 
appears smoother as the number of blocks n increases 
for constant chain length nN as the A-B junctions be- 
tween blocks increase. Of course, the temperature that 
we enter the globule regime (upon lowering the tempera- 
ture) is the same for homopolymer and copolymer chains, 
because in our model the location of this boundary de- 
pends only on the solvent (controlled by temperature in 
our simulations), which is the same for a homopolymer 
chain, as well as for type A and type B beads of a copoly- 
mer chain, but in the case of copolymer chains it only 
appears smoother. 

Further on our simulation method, the positions ri{t) 
of the effective monomers with label i evolve in time t 
according to Newton's equation of motion, amended by 
the Langevin thermostat 
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where Ui{rj{t)) is the total potential acting on the I'th 
bead due to its interactions with the other beads at sites 
rj{t), 7 is the friction coefficient, and Ti{t) the associated 
random force. The latter is related to 7 by the standard 
fluctuation-dissipation relation 
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Following previous work [21|-|30|, we choose 7 
the MD time unit tlj = (wcrlj/egj)^^^ — 1 being also 
unity for our choice of units. Eq uation ^ is integrated 
using the leap-frog algorithm J3l|] , with a time step At = 
0.006Tt,7 and utihzing the GROMACS package ^. 

Simulation of polymer chains under poor solvent con- 
ditions is notoriously difficult as relaxation times become 
exceedingly high. Therefore, we describe here briefly our 
simulation procedure. After equilibration at a "high" 
temperature (in the coil regime), we collect a number 
of independent samples (typically 500), which we use as 
initial configurations for slow cooling runs, as is done in 
previous studies [27h29| |. For longer chains, temperatures 
higher than T = 3.0 were used in order to facilitate the 
procedure of obtaining initial independent samples. We 
note here that the solvent is taken into account in our 
model only implicitly by tuning the temperature |13l - 
Il5| . Then, decrease of the temperature corresponds to 
higher incompatibility of the implicit solvent with the 
monomers. For each cooling history, we lower the tem- 
perature in steps of AT = 0.1 and run the simulations 
for 2 X 10^ MD steps using the final configuration at 
each (higher) temperature as starting configuration for 
the next (lower) temperature. At low enough tempera- 
tures [1J| (typically T < 2.0, at a rather high distance 
from the Q temperature), where dense "clusters" of a 
few neighboring blocks are formed, it is not possible to 
run simple MD simulations long enough to sample the 
phase space adequately. Therefore, using this procedure 
of independent cooling histories is indispensable for ob- 
taining meaningful statistics. This large statistical effort 



prevents the study of exceedingly long chains for this 
simulation model. 

At temperatures close to 0, coil-like structures are 
formed and the differences between multiblock copoly- 
mer chains with different n, N are rather small, expecting 
that at T ^ oo these differences will completely disap- 
pear. At temperatures below T « 2.4 the chains collapse 
and blocks of the same type of monomers form clusters 
with other monomers belonging to different blocks, due to 
the unfavorable interactions between A and B beads |13l - 
[l5j . The static properties of these clusters have been de- 
scribed in detail in previous communication |13j . Here, 
we focus on the static properties of the whole chain (of 
length nN) as well as the static properties of individual 
blocks of N effective monomers. Due to the symmet- 
ric choice of our parameters, the averages of the static 
properties for A and B blocks should not differ. There- 
fore, the calculation for averages of the properties of the 
blocks should be taken over all blocks n, irrespective of 
their type, and their properties will be denoted with the 
subscript "b" , for example the average value of the gy- 
ration radius of individual blocks would read Rgb while 
those for the whole chains without this index, namely as 
Rg. We discuss below the properties that might need 
some more details. 

We follow the definition of standard textbooks |33l - [35| 
to discuss aspects related to the persistence length of the 
polymer chains, where the persistence length is defined 
from the decay of bond orientational correlations along 
the chain. However, one should be careful when defines 
the persistence length with this correlation function as it 
has been discussed elsewhere [3a,|33]. Nevertheless, we 
show in this work that this property is extremely accu- 
rate in indicating the appearance of full phase separation 
of the blocks in symmetric multiblock copolymers. We 
give here a brief definition of this property. Defining the 
bond vectors di in terms of the monomer positions fi as 
di = ri+i — fi, i = 1, ...,nN — 1, the bond orientational 
correlation is defined as 
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FIG. 2: (Color online) Acylindricity of the blocks for 

chains of total lentgh nN ~ 600. The dependence with 

the block length N is shown for different temperatures 

(a) and the dependence with the temperature T for 

different number of blocks n (b). 



Note df = ll, with 4 being the average bond length. For 
this model It w 1. 00 and hence {cos9{Q)) « 1, of course. 
Considering the limit nN — > oo and assuming Gaussian 
chain statistics, one obtains an exponential decay, since 
then {cos6{s)) = (cos0(l))* = exp[ln(cos6'(l))], and thus 
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Below Q temperature this correlation function falls-off 
rapidly at small s depending strongly on temperature T. 
In our study we also consider properties related to the 
shape of multiblock copolymers. We follow the descrip- 
tion of Theodorou and Suter [40| to define the apshericity 



and acylindricity of a mutliblock copolymer chain. First, 
we define the gyration tensor as 
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where s^ = col{xi,yi, Zi) is the position vector of each 
bead, which is considered with respect to the center of 

mass of the beads J2iJi ^* = Oj ^"^^ the overbars denote 
an average over all beads N^. When the gyration tensor 
of the whole chain is considered, then iV^ = nN. For the 
blocks iVj = N, the gyration tensor and the properties 
are calculated for each block separately and then an av- 
erage over all blocks irrespective of their type is taken. 
As mentioned earlier, this is perfectly valid due to the 



symmetry of our model. The gyration tensor is symmet- 
ric with real eigenvalues and a cartesian system that this 
tensor is diagonal can always be found, i.e., 



S = dw5(X2,y^Z2), 



(9) 



where the axes are also chosen in such way that the di- 
agonal elements (eigenvalues of S) X^, Y^, and Z'^ are 
in descending order {X'^ >Y^>Z'^). These eigenvalues 
are called the principal moments of the gyration tensor. 
From the values of the principal moments, one defines 
quantities such as the asphericity &, 



b = X^~l/2{Y^ + Z^). 



(10) 



When the particle distribution is spherically symmetric 
or has a tetrahedral or higher symmetry, then b = 0. The 
acylindricity c 



= y2 _ ^2 



(11) 



is zero when the particle distribution is in sync with a 
cylindrical symmetry. Therefore, the acylindricity and 
asphericity are the relevant quantities that would de- 
scribe geometries that could be relevant in the case of a 
homopolymer or multiblock chain. These quantities are 
taken with respect to S, to the sum of the eigenvalues, 
i.e. the square gyration radius of the chain, which we also 
have calculated independently on our original cartesian 
coordinates in order to check our results. 

Assuming an ellipsoidal shape and based on the calcu- 
lation of the above eigenvalues of the gyration tensor, one 
can define an effective volume expressed by the following 
relation lill 



(12) 



V^ff = 47W^Y[y%, 



where Ai = X'^, A2 = Y'^, and A3 = Y'^, are the eigenval- 
ues of the radius of gyration tensor. Then, the effective 
radius of a sphere with the same volume as this ellipsoid 
(yaff'j is given by the geometrical mean of individual 

This can be compared with 



radii i?^-^/ =\/3nti^A7 

the volume of an effective sphere defined by the gyration 

radius Rg =y'Yyi' ^i- Such an analysis is not in the 
scope of the present study. 

Another property that provides information on the 
microscopic properties of multiblock chains is the num- 
ber of contacts between beads of the same or different 
type. Smaller number of contacts A-B suggests a ten- 
dency towards phase separation. The number of con- 
tacts can be mathematically described by the following 
formula [l3, HI 



„/3 =47r / g„^(Ar)(Ar)''d(Ar), 



(13) 
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FIG. 3: (Color online) The dependence of asphericity 

with the temperature T for different combinations of n, 

N for chains of total length nN — 600 is shown. 



where Ar is the absolute value of the distance between 
two sites of monomers ri , fj in the multiblock copoly- 
mer chain, and ga[j the corresponding radial distribution 
function. Equation (J13p means that a pair of monomers 
[a, /?, (q;,/3 = A,B)] is defined to have a pairwise con- 
tact if their distance is less than r„. For the distance r„ 
we have used the standard Stillinger [43 neighborhood 
criterion for monomers. We followed the standard choice 
r„ = l.StT^j and checked that qualitatively very similar 
results were obtained if one chooses r„ a bit smaller than 
this choice (larger values of r„ are physically hardly sig- 
nificant, since then the particles are too weakly bound, 
due to the rapid fall-off of the LJ potential). There- 
fore, a pair of monomers being an absolute distance less 
than r„ = 1.5 apart define a "contact". Then, the num- 
bers presented in this manuscript will denote the average 
number of neighbors per monomer. 



III. RESULTS 

Firstly, we discuss properties related to the acylin- 
dricity and asphericity of the chains. The acylindric- 
ity defines the deviation of the shape of the chain from 
a cylindrical geometry; correspondingly the asphericity 
expresses the deviation from a spherical shape or tetra- 
hedral or higher symmetry with respect to the total di- 
mensions of the chains. Then Fig. [2^ presents the depen- 
dence of the acylindricity of individual blocks on average 
for multiblock copolymer chains of the same total length 
nN ~ 600 on the block length N . Thus, the number of 
blocks for each chain is n = 600/7V. This choice for the 
total length of the multiblock chain allows for the ob- 
servation of all interesting regimes related to the phase 
behavior of such macromolccules p^3l - (l5j . In Fig. [5^ one 
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FIG. 4: (Color online) (a) Acylindricity of the the 
blocks of a chain with n — 8 blocks as a function of the 
block length N for different temperatures, as indicated. 
The dependence of acylindricity on the number of blocks 

n is shown correspondingly in part (b) for N = 80. 



FIG. 5: (Color online) Asphericity of chains composed 

of n = 10 blocks as a function of the block length N (a), 

and asphericity of chains with block length A^ = 20 as a 

function of the number of blocks n (b). 



can observe very clearly a regime for intermediate values 
of 2 < A^ < 20 where the curves show small deviations 
from one another. The cases of A^ = 2 and A^ = 1 are par- 
ticular. For A^ = 1 (not shown) the acylindricity drops to 
zero. At lower temperatures larger deviations for higher 
values of A^ are seen. Nevertheless from such a plot we 
arc not able to find any hint that relates to the well- 
known behavior of these systems [13l - [la | . Looking at the 
dependence of the above systems with different A^ on the 
temperature T (Fig. [^b) it is clearer that small blocks 
show very small dependence on the temperature due to 
the very small size of the blocks (TV = 600/n), as it also 
happens in the plot of Fig. [2^. The graph of Fig. [5}d in- 
dicates that for small values of A^, i.e. n > 20, the blocks 
rather do not show any dependence on the temperature. 



while for n < 20 (A" > 30) , stronger dependencies on the 
the temperature T are seen. This suggests that stronger 
deviations from a cylidrical symmetry persist for multi- 
block chains of higher A^, as nN remains constant. This 
behavior is also related to the behavior of the size of 
the cluster formations in these systems, which has been 
discussed elsewhere [13|. From the results of Fig. [5] we 
do not see any direct correlation to the phase behavior 
of the multiblock chains by examining the shape of the 
individual blocks. 

On the other hand, the asphericity of the whole multi- 
block chain of total length nN = 600 shows interest- 
ing behavior when a plot with the temperature is at- 
tempted. We remind the reader that the total chain 
length nN = 600 is high enough to consider all the 
different cases of phase behavior of symmetric multi- 
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FIG. 6: (Color online) Asphericity of the blocks of 

niultiblock chains with n = 8 as a function of the block 

length A^ (a), and of chains with block length A^ = 80 

(b) as a function of the number of blocks n. 



block copolymers, which have already been discussed in 
detail [ij, [l^] . Figure [3] then illustrates two different 
regimes. At rather high temperatures (T > 2.4) the dif- 
ferent multibock chains (of the same total length nN) 
exhibit small differences as we are close to & conditions. 
Such small differences are related to the number of prox- 
ima A-B interactions at the junctions that connect the 
consecutive blocks A-B blocks. The values on the j/-axis 
of this plot indicate strong deviation of the chain's shape 
from a spherical one, due to the excluded volume interac- 
tions in combinations with the incompatibility between 
monomers of type A and B at the junctions of A and 
B blocks. At high values of N (small n except for the 
case n — 2 which will be separately discussed below) 
we observe strong differences at low temperatures. The 
case n = 4 reaches at temperature close to T = 2.2 a 



FIG. 7: (Color online) Bond orientational correlation 
for the case n = 10, A^ = 60 (a) corresponding to a full 
phase separated system at temperature below T ~ 2.4 

(here T = 1.5), i.e. all blocks of type A join in one 
cluster while the blocks of type B form another cluster. 

At temperatures T > 2.4 the chain obtains coil-like 
formations. The dependence on the block length A^ for 

chains with n ~ 16 is shown in part (b). For block 

length A^ > 20 the phase separation between blocks A 

and B takes place. These data refer to the temperature 

r==: 1.8. 



plateau-like regime. We can even observe a slight increase 
as temperature decreases and the incompatibility of A 
and B monomers increases. As n increases (TV decreases) 
this plateau-like dependence appears correspondingly at 
lower temperature T. Knowing already the phase behav- 
ior of these systems |13l - [la |. these temperatures corre- 
spond to those that the multiblock copolymer chain sep- 
arates in two clusters with a single A-B interface forming 
between them, where each cluster contains only beads of 



one type with a permanent interface A-B separating the 
two clusters. The transition of multiple A-B interfaces 
to a single A-B interface becomes also noticeable in the 
behavior of asphericity, for high n, i.e. n > 20, where 
this plateau now has disappeared. Therefore, the phase 
behavior of the multiblock chains is nicely reflected to a 
quantity that describes the overall formation of the chain, 
without taking into account the measurement of any mi- 
croscopic quantities which relate to the incompatibility 
between A and B beads. 

For the case n = 2 this transition to the phase sepa- 
rated clusters can not be concluded from the results of 
Fig. |3l because in this case the chain always has only 
these two clusters of A and B beads either in a coil like 
formation at higher temperatures or a globule-like for- 
mation at lower temperatures. In this case of n = 2 we 
rather see a small decrease in the quantity b/S, as the 
chain obtains a globular formation. Also, the values of 
b/S for n = 2 are higher than the other cases indicat- 
ing the stronger presence of A-B interface. In this case a 
smaller number of unfavorable contacts A-B exist (Fig. 4b 
of Ref. [i3|)- For all n shown in Fig. [3] a more spherical 
formation is obtained at lower temperatures. In partic- 
ular this spherical symmetry is strongly present for the 
multiblock chains of small N, where also we do not have 
the appearance of this large A-B interface and therefore 
the curves vary smoothly with the temperature T. The 
above description is valid for all cases where full phase 
separation occurs without restrictions on the length N 
or number of blocks n and reflects totally the full phase 
separation behavior of multiblock copolymers. 

Now we turn our focus on acylindricity and asphericity 
for the blocks and the whole multiblock chains on their 
dependence with the block length N and the number of 
blocks n, so without the fact that results correspond to 
chains of the same total length nN. Figure |4^ shows 
such results for the average acylindricity of individual 
blocks and its dependence on the block length N for a 
multiblock chain containing n — 8 blocks. It has been 
shown that for this number of blocks phase separation 
in two clusters A and B is taking place for high enough 
block lengths (TV > 20) following our simulation proto- 
col. However, the plot of this quantity for individual 
blocks does not reflect the different phase separation be- 
havior for such systems. At temperatures close to the 9, 
i.e. T = 3.0, an almost straight line is seen suggesting 
that chain blocks appear with the same overall formation 
irrespective of their block length N. At lower tempera- 
tures down to T = 1.5 deviations are seen from the above 
behavior and the chains with higher N show a higher de- 
pendence on the temperature T. Larger blocks tend to 
obtain a more cylindrical-like shape as the temperature 
decreases at T = 1.5, something that is seen for all cases. 

However, it becomes clear that even for small N the 
increase of the number of blocks n favors a more cylin- 
drical shape of individual blocks on average. In the case 
that the chain length N is high enough that full phase 
separation can take place {N = 80, Fig. Hb), i-e., one 



DQ 

m 

< 
< 



< 




(a) 



n=2 
n=4 
n=6 
n=8 
n=10 



-©- 

..A- 

■■■■V 



n=12 
n=20 
n=24 
n=30 
n=40 



o- 
-e- 



'n=50 

n=60 

n=100 

n=120 

n=150 



■••■o 



13 ■„ 

e ° ta. 



Fi^^t---g"-a-a-'-fe--a-g-a-&-6 



©Ig-Slg g S 5 n D Q 




(b) 




(c) 



FIG. 8: (Color online) The dependence of the number 
of contacts on temperature for different n for chains of 
total length nN = 600. Contacts A A (a), AB (b), and 

the total number of contacts irrespective of whether 
they are of case A or B (c) are shown. The fraction of A 

and B monomers is the same for cases, as the total 
length of the multiblock chain is kept constant and the 

most symmetrical case is considered, i.e. an even 

number of blocks of the same length alternate along the 

chain. Block lengths fulfilling this condition for 

nN ~ 600 are considered. 



cluster contains the beads of all blocks of type A and 
another cluster beads of all blocks of type B, variations 
due to temperature can become now more apparent. For 
iV = 80 we observe the same tendency that is seen for 
N = 5, but the chains are long enough to observe the 
differences attributed to different temperatures. This 
can also be concluded by examining the n dependence 
of multiblock chains with various block lengths N. The 
graphs of Fig. |4] show that the shape of individual blocks 
is affected by the number and block of lengths, but phe- 
nomena of phase separation can not be identified on an 
individual block-based analysis as a function of block 
length N or number of blocks n. 

Plots showing the overall asphericity of multiblock 
chains as a function oi N oi n (Fig. [5]) neither do reveal 
any entry to the phase separation regime of multiblock 
chains and the behavior is similar to what is seen for the 
average properties of the individual blocks. Aspheric- 
ity shows an independence on n and N for temperatures 
close to Q (i.e. T ~ 3.0), while b/S decays with A^ or 
n obtaining a more globular formation with decreasing 
temperature and increasing n or A^, indicating the entry 
to the regime where the multiblock chain collapses in a 
poor solvent. Showing the corresponding plots of Fig. |4] 
for the asphericity of the blocks (Fig. [6]), we draw the 
same conclusions. However, the asphericity of the blocks 
shows a more pronounced relation to the entry of the 
globule regime. Higher block lengths exhibit a higher de- 
pendence on the change of temperature and decreasing 
the temperature from T = 2.4 to T = 1.5 (Fig. EK) the 
slope of the curve does not change, implying that regard- 
ing asphericity of individual blocks there is a stronger 
dependence on A^, rather than n (Figs.Eb). 

As in the case of Fig. [3] where we have illustrated the 
asphericity of a multiblock chain as a function of temper- 
ature, the phase behavior of multiblock copolymer chains 
can be perfectly reflected to other properties, as for ex- 
ample the orientational correlation function of Eq. ([5|l. 
Then Fig. [7] shows examples of this quantity for different 
cases. For a system not exhibiting full phase separation 
this orientational correlation (not shown here) for high 
values of distance s (defined in number of beads along 
the chain) is usually neglected, since in this range of s 
we rather observe random scattered points indicating a 
loss of correlation. However, as we will see this range of 
s could provide very important information in the case 
of multiblock copolymer chains, when a single interface 
between blocks A and B forms for various n and A^ at low 
temperatures T. For s close to zero (Fig. [7]) one is not 
able to fit any exponents, rather there occurs a very fast 
decay, which is more pronounced at lower temperatures. 
It is remarkable that this initial decay is independent of 
the varied parameters; an obvious interpretation is that 
in the scale of a few subsequent bonds the chain main- 
tains a high local flexibility. 

We focus now on our aim which is to relate the phase 
behavior to the bond orientational correlations, which 
is very well shown in Figs. [7^ and b. Figure [7^ shows 
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FIG. 9: (Color online) The dihedral angle as a function 
of temperature for chains of total length nN — 600 (a). 
The dependence of the same property for chains with 
A = 40 with the number of blocks n (b) , and for chains 
with n = 20 with the block length A^ (c) . D is measured 
in radians. 
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FIG. 10: (Color online) Two representative cases for the 

mean square gyration radius of the blocks for chains 

with n = 20 (a) , and for chains with A^ = 40 (b) for 

different temperatures as indicated. 



results for a multiblock copolymer chain with N = 60. 
At temperatures T < 2.1 and s > 60 we observe peaks 
which are repeated periodically, as s increases. We can 
clearly observe four different peaks, which correspond to 
the (n — 2)/2 consecutive A-B blocks respectively. The 
centers of these peaks correspond exactly to an integer 
multiple of the total length of an A and B block, which 
are joint together one after the other. These peaks ap- 
pear for the range of temperatures that a complete phase 
separation of A and B blocks takes place (a single inter- 
face between A and B blocks is formed). Then the values 
at small s {s < 2N) correspond to correlations of bonds 
within the length of a block A and B. Again, the orien- 
tational correlation function is able to provide very good 
information for intermediate and high s for the phase be- 
havior of multiblock copolymer chains under poor solvent 



conditions. Figure [7|d shows results for the case n = 16 
at T = 1.8 for chains of different block length N. At 
this temperature, systems that exhibit a full phase sepa- 
ration between blocks A and B show symmetrical peaks 
at positions 2/iV, where / takes values 2, 3, . . . , (n — 2)/2. 

The number of contacts of beads docs not also re- 
veal directly any correspondence to the phase behavior 
of multiblock copolymers (Fig. |S]), although these num- 
bers are a direct result of the microscopic interactions 
between monomers. The number of contacts A- A (which 
is equal to the number of contacts B-B due to the symme- 
try of our model) decreases monotonically with increasing 
temperature as beads come apart, but with a lower rate 
for temperatures T > 2.4, as we can see from Fig. [S^. 
We have found [Ij] that this is the temperature that 
the chain enters the collapsed state, i.e. effects due to 
the poor solvent are becoming very important until the 
chain obtains a globular formation at even lower temper- 
atures. For different number of blocks n, this change in 
the slope sets in at different temperatures. Of course, for 
symmetric multiblock copolymer chains of the same to- 
tal length jiN — 600 and increasing n it is expected that 
the number of contacts A-A will become lower. Increas- 
ing the number of blocks and joining them one after the 
other, the number of contacts A-B increases. This is bet- 
ter shown in Fig. [8)d. Counting the number of contacts 
between beads irrespective of their type (Fig. [5}:), we can 
see that the increasing number of blocks n for chains of 
the same total length nN results in smaller amount of 
contacts between beads. This is obviously attributed to 
unfavorable contacts A-B; beads A and B prefer to be 
apart, reducing in this way the number of contacts ntot 
for chains composed of the same amount of beads A and 
B. Again, such plots do not provide a clear information 
on the regime that this phase separation between blocks 
in a multiblock copolymer chain takes place. We rather 
see a change in the slope of the curves upon entering the 
collapsed state, and the results of Fig. [8] suggest that any 
transitions take place gradually. 

As the bond orientational correlation has shown a di- 
rect correlation on the phase behavior of a multiblock 
chain, we have also looked at other properties as a func- 
tion of n, N, and T, as it is for example the dihedral 
angle D which is formed by four consecutive monomers 
and taking the average of this quantity along the chain 
for all nN — 3 tetrads disregarding in this way effects at- 
tributed to different sequences along the chain. As the 
orientational parameters for small correlation distances s 
in terms of monomer count (e.g. s < 4) has shown, the 
chain keeps its local flexibility irrespective of the inter- 
actions A;B, or even the chain architecture |37|. There- 
fore, small differences in D are expected with variation 
of the parameters n, N , and T in the case of multiblock 
copolymers, as it also appears in Fig. ^ Then, Fig. [SJi 
illustrates the behavior of D with temperature T. In- 
crease of the temperature induces a slight increase in D 
for all chains with different number of blocks, that is a 
tiny local stretching of the chains. The increase in the 
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FIG. 11: (Color online) Two representative cases for the mean square gyration radius for chains of n = 10 as a 
function of the block length N (a), and for chains of A^ = 20 as a function of the number of blocks n (b). 



number of blocks n shows that D also increases, but this 
change is rather small due to reasons that have already 
been discussed. In this case all chains have the same to- 
tal length nN = 600. If we choose a rather intermediate 
block length (as for example in Fig. IHb, N = 40), we can 
see the dependence on n. However, in this case and for 
n > 20, D has reached a plateau. For small number of 
n there is a sharper increase of D with n, which is also 
more pronounced at lower temperatures and in the range 
n < 10. But, this is not the whole story, D depends also 
on the block length N (see Fig. |9]:). Here we show an ex- 
ample for n = 20, which corresponds to a point belonging 
to the plateau regime of part (c) of Fig. [3] D provides in- 
teresting information on the local behavior of multiblock 
copolymer chains, with small differences arising by varia- 
tion of the parameters and no effects related to the phase 
behavior of the multiblock copolymers can be found, as 
expected, since D is a particularly local property of the 
chain. 

Figure [TOk shows the dependence on N for a choice 
of intermediate number of blocks n = 20 for different 
temperatures. Differences due to temperature are ex- 
pected to increase as the block length N increases. We 
see that small differences result from the variation of TV 
for different temperatures. Similar behavior we see for 
other values of n. Interestingly the size of individual 
blocks also depends on the number of blocks n (Fig. [TUb) 
at lower temperatures, but at higher temperatures the 
gyration radius of the blocks does not change with n. 
Again, these are results for an intermediate value of block 
length iV = 40, where a full phase separation between A 
and B blocks can take place at low temperatures. One 
could argue that the number of blocks does not play any 
role when the chain has a coil formation, while in the 
collapsed state a n-dependence occurs at low number of 



blocks n. However, a connection of phase separation or 
a partial phase separation between blocks is not seen for 
this property of individual blocks by examining all rele- 
vant cases and possible variations of parameters. Never- 
theless, it shows us that the dependence on n of the indi- 
vidual size of the blocks seen under poor solvent condi- 
tions, is smeared out at high number of blocks n. Similar 
effects are seen for other choices of N. Unfortunately, in- 
creasing the number of blocks n for the interesting range 
of N and apply further exploration of this phenomena is 
rather limited in the present simulations. 

Assuming a rather small number of blocks n, and 
knowing that full phase separation is favored by a small 
number of blocks, we plot the squared average gyration 
radius of the whole multiblock chain (Fig. [TTb). At tem- 
peratures close to 6 one can see a straight line (mind the 
double logarithmic scale), as expected. Deviations from 
this behavior are seen for lower temperatures. The de- 
pendence on the number of blocks n is shown in Fig. Illb . 
In both figures we can rather see the same behavior for 
variation of n and iV, showing that the overall dimensions 
of the multiblock copolymer chain rather do not depend 
very much on the way that the blocks are distributed 
along the chains. Of course, increasing the number of 
blocks favors a more elongated chain at lower tempera- 
tures, due to the higher number of unfavorable interac- 
tions between A and B monomers. Similar conclusions 
are drawn by examining all combinations of n, N , and 
T. 



IV. SUMMARY 

In this manuscript we have presented results for a sin- 
gle multiblock chain and its individual blocks varying 
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parameters such as the block length JV, the number of 
blocks n, and the temperature T. The phase behavior 
of multiblock copolymer chains in a dilute solution has 
been discussed in detail elsewhere [IJj II3 , as well as the 
static properties of the clusters composed of blocks A or 
B [l3|, which are formed due to the microphase separa- 
tion of monomers A and B. In this manuscript we have 
tried to make a connection between the phase behavior 
of these macromolccules and the static properties of the 
chain and of the individual blocks on average. We have 
shown that properties related to the chain as a whole 
can reveal significant information on the phase behavior 
of multiblock copolymers, while study of the individual 
blocks does not provide such information. The presence 
of a single A-B interface between blocks of different type 
play the major role in the behavior of multiblock chains. 
The effect of the presence of this interface is reflected to 
properties such as the bond orientation correlation func- 
tion and the asphericity of the whole multiblock chains. 
In particular, we have shown that correlations for inter- 
mediate and high distances s (measured in number of 
monomers) along the chains for the orientation correla- 
tion function which, is generally not interesting for a ho- 
mopolymer chain, can show very interesting effects due 
to the presence of a well defined A-B interface in the case 
of multiblock copolymers. Furthermore, it is remarkable 
that the study of properties of the individual blocks can 
reveal dependence on parameters such as the number of 
blocks n as the chain gets into the collapsed state. This 
dependence is lost for higher number of blocks. Phase 
separation that is due to the microscopic interactions 
between monomers and the incompatibility of A and B 
monomers can be mainly measured with properties re- 
lated to the chain as a whole, or on the basis of a cluster 
analysis for multiblock copolymer chains, as was done 
previously |13l - [l5| . 

The results presented in this manuscript, in com- 



bination with recent results in the phase behavior of 
multiblock copolymer chains, give a detailed descrip- 
tion for these systems, at least for the most symmetrical 
case, where an even number of blocks of A and B type 
monomers alternate along a linear multiblock chain with 
A^^ = Nb being the lengths of blocks A and B respec- 
tively. It is interesting to sec that the change of com- 
position of such symmetrical chains can lead to different 
scenarios of phase separation. Introducing any kind of 
asymmetry in our systems or the consideration of an ad- 
ditional type of monomers in the structure of multiblock 
chains could result in different and interesting behavior 
for these systems. The use of a standard generic coarse- 
grained model to describe our chains could act as a refer- 
ence system for the study of these complicated systems. 
Moreover, the behavior of multiblock copolymers can be 
parallelized to that of various biological macromolccules 
which are formed by periodically repeated chemical units 
("monomers") along their chain. A detailed discussion 
on the properties of multiblock copolymer chains using a 
simple model for a variety of parameters n, N^ and T is 
useful to understand the behavior of these chains, before 
any other effects, due to, e.g. charges, more complicated 
interactions, architecture, etc., come into play. We en- 
visage that the present contribution could also stimulate 
further work on these topics. 
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